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For transport processes in geometrically restricted domains, the mean first-passage time (MFPT) 
admits a general scaling dependence on space parameters for diffusion, anomalous diffusion, and 
diffusion in disordered or fractal media. For transport in self-similar fractal structures, we obtain 
a new expression for the source-target distance dependence of the MFPT that exhibits both the 
leading power law behavior, depending on the Hausdorff and spectral dimension of the fractal, as 
well as small log periodic oscillations that are a clear and definitive signal of the underlying fractal 
structure. We also present refined numerical results for the Sierpinski gasket that confirm this 
oscillatory behavior. 
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The time taken for a random walker to reach for the 
first time a given target site from a given source site, 
is a particularly interesting physical observable that en- 
codes information about a wide variety of physical and 
mathematical transport processes, ranging from chemical 
transport to propagation through networks [1]. Recent 
progress has led to precise quantitative expressions that 
show how the mean first passage time (MFPT) is af- 
fected by the geometry [e.g., volume, shape, dimension] 
of the region in which the random walk occurs, even for 
fractals [2, 3]. These results have been confirmed by nu- 
merical simulations and have motivated various exper- 
imental studies [4]. Here we argue that for such pro- 
cesses there are observable finer details encoded in the 
MFPT when the transport occurs on highly symmetric 
fractals: those possessing exact self-similar scaling, such 
as the Sierpinksi gasket. Technically, this result requires 
information concerning the off-diagonal Green's function 
G{yt\?s) between the source and target points, and we 
show that for highly symmetric fractals, simple nearest 
neighbor random walks gencrically lead to spatial log- 
periodic oscillations in first passage observables such as 
the MFPT. 

We recall briefly the essential results from [2, 3]. Con- 
sider a random walker moving in a bounded region of N 
sites. In [2, 3] it has been shown that the MFPT to a 
target site starting from rs satisfies in the large N 
limit: 

(T)/iV~G(r T |r T )-G(r T |r s ). (1) 

Note that the r.h.s. is independent of the confinement 
and involves the infinite space Green's function G(r|r') 
defined by 

r-OO 

G(r|rO= / P(r,tW)dt, (2) 
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where the propagator P(r, £|r') gives the probability den- 
sity that a random walker evolving in infinite space is 
at site r at time t, starting from site r' at time 0. 
This probability density satisfies the diffusion equation 
d t P(r,t\r') = A r P(r,i|r'). 

We show in this paper that for certain highly symmet- 
ric fractal structures, the Green's function when averaged 
( 7TT ) over source and target points, has the general scal- 
ing form: 

where r — \tt — i*s| is the source-target distance, and I is 
a decimation constant specific to the fractal (see Table 1 
in [5] for some examples). The function Q(x + 1) =Q(x) 
is a periodic function of period unity. The fractal (Haus- 
dorff) dimension df and the walk dimension d w [6, 7] are 
related to the spectral dimension as d s = 2df/d w . The 
walk dimension d w accounts for the anomalous space- 
time dispersion on fractals, (r 2 (t)) <~ t 2 / dvj , generalizing 
the usual expression, (r 2 (t)) <~ t, for Euclidean manifolds 
[24]. The scaling behavior of G(tt\ts) in (3) has in- 
teresting consequences for a number of relevant physical 
quantities. Adapting the arguments of [2, 3], we deduce 

r A-r d -- d fg(2TT^j) , d w <d f 

(T)/iV~ I A + \nrg(2n^) , d w = d f (4) 
[ r d w -d f g ( 2?r W) ; dw>df 

where A is independent of both TV and r. Previous work 
[2, 3] took Q = c, a constant, and found excellent agree- 
ment with numerical results for a wide variety of diffu- 
sion processes: sub-diffusive (d w > df), super-diffusive 
(d w < df), and critical (d w — df). 

However, closer inspection of the numerical results in 
[2, 3] indicates that for certain highly symmetric fractals 
there are small oscillatory deviations from the leading 
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forms from (4) with Q = c. We argue here that these 
oscillatory deviations are not numerical artifacts, but are 
in fact real physical results that are strong indications of 
the underlying scaling and fractal structure of the diffu- 
sion process. We find an explicit approximate form for 
the function Q which quantitatively explains these small 
oscillatory effects for the Sicrpinski gasket. 

It is clear from (1) that we need information about 
the spatially off-diagonal propagator. The approxima- 
tion of taking Q — c coincides with assuming a scaling 

form P(r,i|r') = IT (^Q-) /t d f/ d ™ [8]. But we can im- 
prove on this estimate. Upper and lower bounds for the 
propagator P(r,t\r') have been proven by Barlow and 
Perkins [9] for a wide class of diffusion processes: 

P t ( Cl ,c 2 ,r,r') <P(r,i|r') < P t (c 3 , c 4 , r, r') (5) 

with 



T t (ci,c v r,r') = exp 



-a 



(6) 



for real and positive constants c\, . . . , C4. The bounds 
(5) translate into corresponding upper and lower bounds 
for the Green's function (here for df 7^ d w ): [25] 
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r — r 



l\dt—d u 



< G{r\r') < 



r — r 



l\df—d u 



(7) 



Unfortunately, little is known rigorously about the con- 
stants c\,---,cq, or about how P(r,t\r') or G(r|r') be- 
have between the respective bounds (5) and (7). 

In contrast to these spatially off-diagonal quantities 
P(r, t\r') and G(r|r'), much more is known about the spa- 
tially diagonal propagator, P(r,i|r) = (r|e tA |r), (which 
when integrated over space gives the heat kernel trace): 



P(r,t\v) = 



1 /2nlnt 
t d sl d ™ \d w \nl 



(8) 
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The function F(x) = F(x + 1) is periodic of period unity, 
which often is well approximated by its first harmonic 
[5], and I is the same decimation constant appearing in 
(3). The log-periodic oscillations in the time variable t 
are ubiquitous and already well recognized not only for 
diffusion on fractals, but also in other complex systems 
having a discrete (lacunar) scaling symmetry [10]. The 
origin of these log-periodic oscillations in t can be traced 
directly to spectral properties of the Laplacian A. To see 
this, write P(r,£|r) in terms of the associated zeta func- 
tion Ca(s) defined by the inverse Mellin-Laplace trans- 
form [26]: 



tr r (e* A ) = ^-Jj- s T( S )(: A {s)d S . 



(9) 



For symmetric fractals [5], Ca(s) has a tower of poles at 
s n = df /d w +2Tvin/(d w In/) in the complex s plane, which 
lead to the log-periodic oscillatory form in (8). Physi- 
cally, these complex poles are a consequence of the expo- 
nential growth of the degeneracies and eigenvalues for the 
Laplacian on a fractal [11]. To have a better understand- 
ing of how this works, it is interesting to go back to the 
familiar case of a regular rf-dimensional Euclidean mani- 
fold [12, 13] where the Laplacian eigenvalues have instead 
a polynomial growth, A„ ~ n 2 , while the degeneracy fac- 



tor scales like deg n 



~ n 



d-l 



for a symmetric manifold, 



such as a d-dimensional hypersphere [14]. Thus the zeta 
function goes like £a(s) ~ Tl in nd ~ 1 l n2s ~ (,r{2s — d+1), 
(£r is the Riemann zeta function) which has the familiar 
real pole at s = d/2. This leads, via the Mellin transform, 
to the well known behavior of the diagonal propagator, 
P(r,i|r) - J2n nd ~ le ~ n2t ~ lA d/2 , at short times. 

On the other hand, for a symmetric fractal, with ex- 
ponentially growing degeneracies and eigenvalues, 



deg„ 



(10) 



we find C(s) - £„ a n /b ns - 1/(1 - a/b s ), which has a 
tower of complex poles with real part s = lna/ln&, and 
vertical spacing 27r/ln6, which identifies In b — d w In/, 
and lna = df In/, so that lna/ln6 = df/d w . The 
complex poles produce the log-periodic oscillations in 
(8). For such a fractal, P(r,t|r) - E„a"e" h "* ~ 
P(^lni) /t lna / lnb , where F is a periodic function of 
ln£, as in (8). For mathematical discussions of oscilla- 
tions in heat kernel estimates see [15, 16]. In particular, 
there is a class of fractals where oscillations are related 
not to the high degeneracies of eigenvalues but rather to 
large gaps in the spectrum. This topic is a subject of 
active current research [17, 18, and references therein]. 

Since space and time are coupled scaling variables re- 
lated through the walk dimension d w , we might expect 
to observe an analogous type of log-periodic oscillations 
for the stationary Green's function G(r|r'), as a function 
of the spatial distance between the points. To see this, 
consider an eigenfunction expansion: 



G(r|r')=£deg, 



0n(r / )0n(r) 

An 



(11) 



On a regular Euclidean manifold, when averaged (~) 
over the points r and r' we find </>„(r') <fi n (r) ~ f(r/L n ), 
where r denotes the distance |r — r'|, and L n depends on 
n. Consistency with scaling and conservation of probabil- 
ity determine L n ~ 1/n, using the polynomial growth of 
degeneracies and eigenvalues, as discussed above for the 
diagonal propagator. We recognize 1/L n = ^/Xn as the 
momentum, reflecting the usual quadratic dispersion re- 
lation, k n ~ \fKi- On a fully symmetric (or unbounded) 
manifold such as a c?-dimensional hypersphere, this be- 
havior holds without averaging. If F(s) is the Mellin- 
Laplace transform of f(r), then /(r) = ^ J c dsF(s)/r s , 
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and 



G(r,r') 



,d-3 



2tti 
ds F(s) 



F(s) 



C \ nr ) 



Cr(s - d + 3) 



1 



,,d-2 



(12) 



iq 2wi 

which is the familiar form of the Green's function, for 
d =/= 2, and is seen to result from the pole of the Riemann 
zeta function at s = d — 2. 

To generalize this result to fractals, we note the em- 
pirical numerical result that averaging over target and 
source points on fractals also leads to an averaged form 
</>* (r) </> n (r') ~ f(r/L n ), for some L n . Consistency with 
scaling and the conservation of probability determines 
L n = b~ n l dw = aT n l d s = l/l n , now using the exponen- 
tial growth of degeneracies and eigenvalues, as discussed 
above for the diagonal propagator. Thus we find a scal- 
ing relation that involves the anomalous walk dimension 
d w , and (for the non-critical case, df ^ d Wl i.e. d s ^2): 



G(v T \v s ) 



Ids P{S) 
b n 2ni J r ) 

F{s) 1 



1 

2wi 



ds ■ 



1 



(13) 



b i+«/d„ 



which leads immediately to the scaling in (3). 

When comparing (3) with the time-dependence of the 
diagonal propagator in (8) , we notice that the arguments 
of the periodic functions differ by a factor l/d w , which 
is consistent with the anomalous scaling of distance and 
time for diffusion on a fractal. Thus, we argue that the 
scaling in (3), and the associated log periodic oscillations, 
have the same physical origin as their temporal counter- 
part in the diagonal propagator P(r, t\r) in (8), coming 
from the exponential behavior of the degeneracies deg„ 
and eigenvalues A„. 

The result (13) can be viewed as a refinement of the 
Green's function bounds in (7), in a similar sense to the 
refinement (8) of the diagonal propagator bounds in the t 
variable in (5, 6). The expression of G(r|r') in (13) is con- 
sistent with (but more explicit than) the Barlow-Perkins 
bounds (5). It is also consistent with rigorous results 
in [19-21] for the small t behavior of t d f/ d ™ lnP(r,i|r'), 
which behaves also as some log-periodic function of r. 

The expression (3) for the averaged stationary Green's 
function constitutes the main result of this Letter. Wc 
now present a numerical test of its consequences for the 
mean first passage time (MFPT) (T). A number of re- 
cent works [2, 3] have derived precise expressions for the 
MFPT in the limit of large domains. Motivated by the 
case of the diagonal propagator (or heat kernel) studied 
in [5], we propose to approximate the log periodic func- 
tion Q in (3) by its first harmonic: 



1 + 62 cos 



2n In'/ 
In/ 



(14) 



Inserting this expression for Q into (4) we obtain new 
expressions for the MFPT (T), which generalize those 
found in [2, 3]. In Figures 1 and 2, we compare this 




FIG. 1: Plot of the spatial dependence of the MFPT ob- 
servable (T) /N, as a function of the source-target distance 
r, for the Sierpinski gasket. The dotted points are numeri- 
cal data (gaskets of generation 3,4,5,6 respectively in yellow, 
blue, green, red), the dashed black line is the leading behav- 
ior r dw ~ dh , as studied in [2, 3]. The blue solid line is a more 
refined fit from (14) with 62 = 0.05, and cj> = n/2. 
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FIG. 2: Plot of the corrections to the leading r dw ~ dh behavior 
of the spatial dependence of the MFPT (T) /N, as a function 
of the source-target distance r, for the Sierpinski gasket. Cir- 
cles stand for numerical simulations (gaskets of generation 
3,4,5,6 respectively in black, red, green, blue), and the blue 
solid line is a fit from (14). This figure clearly shows the 
log-periodic oscillations with the variable 27rlnr/ln/. 

new expression for the MFPT with numerical data, for 
the Sierpinksi gasket, which has df > d w . First, in Fig- 
ure 1 we use the numerical results of [2, 3], and compare 
with the leading power-law behavior with Q = 1, as well 
as with our refined log periodic oscillation form in (14). 
We have used a simple fit with a single periodic function, 
with period In / in the variable lnr. The agreement with 
this refined form is very good. In Figure 2, we present 
the result of a new and more detailed numerical analysis 
for the corrections to the leading behavior, plotted as a 
function of the logarithmic separation distance lnr. This 
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plot exhibits clear log-periodic oscillations of correspond- 
ing period hi I = In 2 for the Sierpinski gasket. 

As a further test of our analysis, one can note the ab- 
sence of oscillations in the case of fractal trees (e.g. the T 
graph studied in [3]). This is related to the fact, proven 
in the mathematics literature, that on trees the Green's 
function is essentially linear with respect to the effective 
resistance distance, in which case no spatial log-periodic 
oscillations are expected. The Green's function and effec- 
tive resistance on trees were studied in [22] , and on more 
general fractals in [23] . The absence of oscillations was 
also checked numerically in the case of disordered frac- 
tals on the example of critical percolation clusters (see 
Figure 3 of [3]), which strongly suggests that their exis- 
tence requires an exact decimation symmetry, found in 
deterministic fractals such as the the Sierpinski gasket. 

In conclusion, we have argued that a refinement of pre- 
viously obtained bounds for the stationary Green's func- 
tion G(r|r') associated to the Laplacian on symmetric 
fractals can be found that is given by expression (13). 
In addition to a leading power law dependence already 
found in [2, 3], we have shown that G(r|r') exhibits also 
small log-periodic oscillations in the spatial variable. Our 
argument traces the physical origin of this refined func- 
tional form to the exponential growth of degeneracies and 
eigenvalues for the Laplacian on a symmetric fractal, and 
is consistent with the anomalous dispersion for diffusion 
on a fractal. The resulting spatial form parallels the tem- 
poral behavior of the diagonal heat kernel, which is well- 
studied both physically and mathematically. While much 
less is known rigorously for the spatial dependence, our 
conjectured form has an immediate consequence for the 
mean first-passage time (T), which is accessible to nu- 
merical and even experimental investigation. We have 
shown that our refined expression provides a good quan- 
titative fit which also demonstrates unambiguously the 
existence of log-periodic oscillations. It would be very 
interesting to investigate the existence of similar behav- 
iors for diffusion on other classes of complex scale-free 
networks, beyond the leading order results in [2, 3]. 
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